Multi-Component Bell Inequality and its Violation for Continuous Variable Systems 
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Multi-component correlation functions are developed by utilizing d-outcome measurements. 
Based on the multi-component correlation functions, we propose a Bell inequality for bipartite 
d-dimensional systems. Violation of the Bell inequality for continuous variable (CV) systems is 
investigated. The violation of the original Einstein-Podolsky-Rosen state can exceed the Cirel'son 
bound, the maximal violation is 2.96981. For finite value of squeezing parameter, violation strength 
of CV states increases with dimension d. Numerical results show that the violation strength of CV 
states with finite squeezing parameter is stronger than that of original EPR state. 
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I. INTRODUCTION 

In their famous paper of 1935 pj, Einstein, Podolsky 
and Rosen (EPR) questioned the completeness of quan- 
tum mechanics, based on a gedanken experiment involv- 
ing the position and momentum of two entangled par- 
ticles. Einstein believed that there must be elements of 
reality that quantum mechanics ignores. It is argued that 
the incomplete description could be avoided by postulat- 
ing the presence of hidden variables that permit deter- 
ministic predictions for microscopic events. Furthermore, 
hidden variables could eliminate concerns for nonlocality. 
For a long time, EPR argument remained a philosophi- 
cal debacle on the foundation of quantum mechanics. In 
1964, John Bell made an important step forward in this 
direction by considering a version based on the entan- 
glement of spin-1/2 particles introduced by Bohm. Bell 
2] showed that the assumption of local realism had ex- 
perimental consequences, and was not simply an appeal- 
ing world view. In particular, local realism implies con- 
straints on the statistics of two or more physically sepa- 
rated systems. These constraints, called Bell inequalities, 
can be violated by the statistical predictions of quantum 
mechanics. 

Although most of the concepts in quantum informa- 
tion theory were initially developed for quantum systems 
with discrete variables, many quantum information pro- 
tocols of continuous variables have also been proposed 
In recent years, quantum nonlocality for position- 
momentum variables associated with original EPR states 
has attracted much attention. The original EPR state is 
a common eigenstate of the relative position i\ — Xi and 
the total linear momentum pi + p2 and can be expressed 
as a 5-function: 



&(xi,x 2 ) = 



e (2Tri/h)(x 1 -x 2 +xo)p ( ip_ 



(1) 



It is important to choose the appropriate type of observ- 
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ables for testing nonlocality for a given state. Bell has 
presented a local realistic model for position and momen- 
tum measurements based on spin-1/2 particles. He has 
also argued that the original EPR state would not exhibit 
nonlocal effects since the Wigner function representa- 
tion of the original EPR state is positive everywhere and 
therefore admits a local hidden variable model. Recently, 
Banaszek and Wodkiewicz Q invoked the notion of par- 
ity as the measurement operator and interpreted the 
Wigner function as a correlation function for these par- 
ity measurements. They then showed that the EPR state 
and the two-mode squeezed vacuum state do not have a 
local realistic description in the sense that they violate 
Bell inequalities such as the Clauser and Home inequal- 
ity |5| and the Clauser-Horne-Shimony-Holt (CHSH)R| 
inequality. The starting point of the demonstration in [4| 
is that the two-mode Wigner function can be interpreted 
as a correlation function for the joint measurement of the 
parity operator. In the limit r — > oo, when the original 
EPR state is recovered, a significant violation of Bell in- 
equality takes place, however, the violation is not very 
strong. To avoid the unsatisfactory feature, Chen et al. 
introduced "pseudospin" operators based on parity, 
due to the fact that the degree of quantum nonlocality 
that we can uncover crucially depends not only on the 
given quantum state but also on the "Bell operator" [j| . 
From reference , the violation of CHSH inequality for 
the original EPR state can reach the Cirel'son bound 
2V§. 



In this paper, we propose a Bell inequality, which is 
based on multi-component correlation functions, for bi- 
partite systems by utilizing d-outcome measurements. 
We then investigate violation of the inequality for contin- 
uous variable systems. Due to the considered d-outcome 
measurements, violation of original EPR state can ex- 
ceed Cirel'son's bound, the maximal vioaltion is 2.96981. 
The CV case with finite value of squeezing parameter is 
also studied. The violation strength of CV states with fi- 
nite squeezing parameter is stronger than that of original 
EPR state. 
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II. BELL INEQUALITY FOR 
MULTI-COMPONENT CORRELATION 
FUNCTIONS 

We consider a Bell-type scenario: two space-separated 
observers, denoted by Alice and Bob, measure two 



different local observables of d outcomes, labelled by 
0, 1, N(= d — 1). We denote Xi the observable mea- 
sured by party X and x% the outcome with X = A,B(x = 
a,b). If the observers decide to measure Ai,B 2 , the re- 
sult is (0, 4) with probability P(ai = 0, b 2 = 4). Then let 
us introduce d A-dimcnsional unit vectors 



v = (1,0, 0,0,--- ,0,0) 

Vl : 



v 2 



1 VN 2 - 1 



N' N 



0,0,- ■■ ,0,0 



1 



1 (N + 1)N N-2 



(N+l)N 



N 7 NVN(N-l)' N V(N-l)(N-2) 



,0,--- ,0,0 



VjV-1 



vat 



1 1 {N + 1)N 1 {N+1)N 
N'~N\ N(N- l)'~iV V (AT - 1)(N- 2)' 



1 (N + 1)N 1 



(N+1)N 



N ' NVN(N-I)' Ny (N — 1)(N — 2)' 



1 

TV 



3-2 ' N 



1_ (N + l)N 
' N\ 3-2 ' 



/(AM 


- l)iV 






1 (N + 


l)iV N 







(2) 



These d vectors satisfy following properties: 



AT 

w E 

j=0 



v,=0 



(u) 



v j' v fc = -^ (j^k) (3) 



For = 2, it is just two valued variables (i.e., v a — 1, V\ — 
— 1) obtained from a measurement. If the measure result 
of Alice is to, and Bob's result is n (where to and n are 
less than N), we then associate a vector v m+ „ for the 
correlation between Alice and Bob [v m+n understood as 
v t , where t — (to + re), mod d]. Based on which, we now 
construct multi-component correlation functions: 



Qij = y~] v m + n P(ai = m, bj = re) 

m.n 

d 

= ^v t P{m + n = t) (4) 
t=o 

where P(a, = m, 6j = n) is the joint probability of 
obtain outcome to and 6, obtain outcome re 



and = 

><5if _1) ) 5 Qy } represents the fc-th 
component of the vector correlation function Qij. 



We now define a Bell quantity for the multi-component 
correlation functions, 



B d 



B (o) 



+ ■■■ + 



(JV+1)JV 



'(iV- l)(iV- 2) 
(A + 1)AT 



B (2) 



2 • 1 



(N + 1)N 



B (N-1) 



N-l 



= E 



k=0 



'(N + l-k)(N-k) 
(N + 1)N 



B (k) 



(5) 



where 



B (0) 



Q\2 



)(0) 
'22 



s<*> = q^-q^-q^ + q^, (mo) (6) 

Any local realistic description of the previous Gedankcn 
experiment imposes the following inequality: 



2(S 2d + (l-6 2d )- 



1 



1 



<B d <2 



(7) 



Obviously, this inequality reduces to the usual CHSH in- 
equality for d = 2. 

The quantum prediction for the joint probability reads 

P QM (ai = to, bj =n) = (t/j\P{a t = to) ® P(6 3 = n)|V)(8) 

where i,j = 1,2; m, re = 0,...,N, P(a, = to) = 
U A \m){m\llA is the projector of Alice for the i-th mea- 
surement and similar definition for P(bj = n). 
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It is well known that the two-mode squeezed vacuum 
state can be generated in the nondegenerate optical para- 
metric amplifier (NOPA)jjj is 



|NOPA) = e r ^ a i- aiO2 )|00) = J2 



(tanhr) r 
coshr 



|nn), (9) 



where r > is the squeezing parameter and \nn) = \n) 1 ( 



\n) 2 = ^{a\) n (al) n |00). The NOPA states |NOPA) can 
also be written as Q: 

|NOPA> = V 1 - tanh 2 r / dq J dq'g(q, q';tznhr)\qq'), 



q 2 +q' 2 —2qq'x 



2(l-x 2 ) 



and 



where g(q,q';x) = -=i — = exp 
V'l 1 - 1 ) 

\qq') = \q) 1 (g) \q') 2 , with \q) being the usual eigen- 
state of the position operator. Since lim x ^i g (q, q';x) — 
S(q-q'), one has linv-.oo J dq f dq'g(q,q';t&nhr)\qq') = 
J dq\qq) = |EPR), which is just the original EPR state. 
Thus, in the infinite squeezing limit, |NOPA) | r _ voo be- 
comes the original, normalized EPR state. Following 
Brukncr et al. .1C| . we can map the two-mode squeezed 
state onto a (i-dimensional pure state: 



sechr 



yl — tanh 2d 



^(tanhr) 1 



(10) 



r n=0 



If the measurement result of Alice is m photons, and 
Bob's result is n photons, we then ascribe a vector 
v m+n for the correlation between Alice and Bob. And 



P(c 



n) is the joint probability of ai obtain 



Local realistic description imposes —4 < B < 2. Nu- 
merical results show that B d=3 (r = 1.4068) ~ 2.90638; 
B d=3 (r -> oo) = 4/(6a/3 - 9) ~ 2.87293. For B d=3 {r -> 
oo), the four optimal two-component quantum correla- 
tions read: 



m photons and bj obtain n photons. More precisely, for 
the two-mode squeezed state one obtains following joint 
probability 

P QM ( ai = m, bj =n) = (ip d \P(ai = m) <8> P{b 3 = n)|V#l) 



III. SOME EXAMPLES 

For d = 3, we have three outcomes vo = (1,0), 
vi = (-l/2,V3/2), v 2 = (-1/2,-^3/2). Accordingly 
the NOPA state is divided into three groups, namely, 

1 00 / 

INOPA) = V tanh 3rl r|3n)|3n) 

cosh r V 

n=0 v 

+ tanh 3n+1 r|3n+ l)|3n+ 1) 

+ tanh 3n+2 r|3n + 2)|3n + 2)j (12) 

The two-component correlation functions depend on 
quantum version of joint probabilities 

P QM ( ai = m,bj = n) 

= (NOPA|P(a 4 = m) <g> P(bj = n)|NOPA) (13) 



Qn — Q22 — QI2 — ( 

1 2, 



3' 3 ; 



2\/3 + l 2-V3, 



V5 
3 



m? = -v. 



(14) 



One thing worth to note that one can treat those three 
two-dimensional vectors in terms of complex numbers, 
namely, vq = 1, v% = uj and vi = uj 2 for simplicity, where 
uj = exp[i27r/3]. Now the Bell inequality becomes 

- 4 < Re[Qu + Qia - Q 2X + Q22] 

+ ^Im[Q 11 -Q 12 -Q 21 +Q 22 ] <2. (15) 

In this sense, we can check the generalized "parity" op- 
erator (cu) n other than usual parity operator (—1)™, the 
two-component correlation function also reads 

Qij = (N0PA|Z<4®Z4 H" a+flb U A ®U b |NOPA) (16) 

where Ua,b are generally U (3) transformations and we 
can sufficiently take them as the products of three spin- 
coherent operators: 

U A = e ^-^°- eb*+-&v- e iJ + -eJ- 

where £j = ^-e^ lVi [actually, the phases <pj can be set 
to be zero, since they do not affect the maximal viola- 
tion. Hence Ua,b is a 50(3) rotation]. I±, U± and V± 
are pseudo-su(3)-spin which can be realized by the Fock 
states as 

00 00 
/+ = |3n)(3n+ 1|, L = £ |3n + l)(3n|, 



n=0 
1 



n=0 



I, = |(|3n)(3n|-|3n + l)(3n+l|), 

OO OO 

f7 + =2|3n + l)(3n + 2|, = £ |3n + 2)(3n + 1|, 

n=0 n=0 

C> 2 = -(|3n + l)(3n + l| - |3n + 2)(3n + 2|) 

OO OC 

y+ = ^|3n)(3n + 2|, VL = £ |3n + 2)(3n|, 



14 = — (|3r?,> <3r?,| - |3n + 2) (3n + 2|). 

Operators {/±,/ z }, {£^±,£4} and {V±,V^} form three 
SU(2) groups, respectively, and {I±,U±,V±, I z , (U z + 
Vz)/v3} forms a SU(3) group. 

We can similiarly get B d (r = finite value) and B d (r — > 
00) with different d. We list them partly in Table [I] 
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Obviously, the degree of the violation increases with di- 
mension d, and the violation strength of CV states with 
finite squeezing parameter is stronger than that of origi- 
nal EPR state. When squeezing parameter goes to infin- 



ity, i.e., original EPR state gotten, we can find the bound 
of the violation of multi-component correlation-function 
Bell inequality [jj . 



[d/2M 2k „ 1 



HmS d (r-oc) = lim4d g (1 - — + , ~ _ , + W 



■%T*[n TT7^-T, ^TTTv?] — 2.96981 (17) 

7r 2 £? (k + 1/4) 2 (fc + 3/4) 2J v ; 



(Ba) 


d= 5 


d= 10 


d= 15 


d= 20 


d = 25 


EPR 


2.91055 


2.9398 


2.94973 


2.95473 


2.9577 


NOPA 
with r 


2.9886 
(1.44614) 


3.03842 
(1.72082) 


3.06836 
(1.8366) 


3.08273 
(1.96562) 


3.08932 
(2.07377) 



equality for bipartite d-dimcnsional systems is developed 
accordingly. We then investigate violation of such Bell 



TABLE I: Violation of multi-component Bell inequality for 
|NOPA) and |EPR) states with different d. 



It is interesting to note that for Bd(r — ► oo), the four 
optimal multi-component quantum correlations share the 

same module: 



2d— 1 



When d tends to infinity, 



\Qij \ ~ V 3d 

\Qij\ = \f2~jz. However, we do not have an analytical 
way to find a bound for violation with finite squeezing 
parameter. For this case, what we do is draw a graph to 
see the variation of Bd(r = finite value) with increasing 
dimension, see Fig^ We calculate the maximal quantum 
violation for CV states with different d. The more the 
value of dimension, the more difficult to find a maximal 
violation. Hence the violation strength points we get are 
for d < 330. With these values, it is easy to see that the 
violation increases from slowly to slowly with increasing 
of d. Which means that there exists a limit for quantum 
violation when d goes to infinity. Until now, we do not 
have an exact value of the limit. We use a software, which 
can give experimence expression given enough points, to 
find a expression that describes the curve in FigQ] 

B = 3.12885 - 1.06535/d + 2.13122/d 2 - 2.19262e~f 18) 

When d — ► oo, quantum violation (B) goes to 3.12885. 
Hence, such value can be thought as an approximate vio- 
lation limit for CV states with finite squeezing parameter. 



IV. CONCLUSION 

In summary, we construct multi-component correlation 
functions based on d-outcomc measurements. A Bell in- 




150 200 

Dimension d 



FIG. 1: Violation of multi-component Bell inequality for CV 
states with finite squeezing parameter for different d. 



inequality for continuous variables case. The degree of 
the violation increases with dimension d, and the limit 
of the violation for the original EPR state is found to 
be 2.96981, which exceeds the Cirel'son bound. The rea- 
son for this is due to the fact that we consider d(> 2)- 
outcome measurements. Numerical results show that the 
violation strength of CV states with finite squeezing pa- 
rameter is stronger than that of original EPR state. 

This work is supported by NUS academic research 
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nancial support from Singapore Millennium Foundation 
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Note added: While completing this work we learn a 
similar result obtained in Ref. which based on the 



CGLMP inequality [U 
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